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We introduce phase operators associated with the algebra su(3), which is the appropriate tool to 
describe three-level systems. The rather unusual properties of this phase are caused by the small 
dimension of the system and are explored in detail. When a three-level atom interacts with a 
. . . quantum field in a cavity, a polynomial deformation of this algebra emerges in a natural way. We 

04 , also introduce a polar decomposition of the atom-field relative amplitudes that leads to a Hermitian 

. relative-phase operator, whose eigenstates correctly describe the corresponding phase properties. 

' We claim that this is the natural variable to deal with quantum interference effects in atom-field 

(N ■ interactions. We find the probability distribution for this variable and study its time evolution in 

^ . some special cases. 

(D ■ 

Q ' PACS numbers: 42.50.Ct, 42.50.Dv, 42.50.Hz, 42. 50. Fx 

(N : 

I. INTRODUCTION 

>: . . . 

: Atomic coherence is essential for a proper understanding of many important effects appearing in the response of an 
T— I ' atomic system to strong laser radiation. Perhaps the MoUow spectrum for a strongly driven two-level atom has 
\ been the first milestone in quantum optics where atomic coherence plays a major role. 

Once the population dynamics in two-level systems was understood, the interest in multilevel (and, especially, 
three-level) atoms became soon apparent. A significant amount of research has thus been devoted to their study 
, and many effects based on quantum coherence such as spectroscopic dark states |^, electron shelving Q, narrow 
spectral lines [Q, pulse matching and anti- intuitive excitation have been highlighted. These nonclassical 
features have an enormous variety of interesting and nontrivial consequences, including electromagnetically induced 
transparency [Q, lasing without inversion state-selective molecular excitation and demonstrations of slow 
light [|lO| and fast light to mention only a few examples. 
^ ■ For a three-level atom coupled to one or two quantized cavity field modes both the atomic dynamics and the 
^ ' statistical properties of the reemitted field have been discussed in detail Jl2| showing interesting phenomena, the most 
outstanding being, perhaps, the existence of collapses and revivals. The phase properties of this model have also been 
studied, mainly u sing the Pcgg-Barnctt formalism for the field, and the connection between revivals and phase has 
been put forward fl^ . 

' In spite of these achievements, even a cursory look at any classical model of atom-field interaction immediately 
H \ reveals that the natural way of understanding the resonant behavior is in terms of the relative phase between the field 
and the atomic dipole. It is also clear that this relative phase is basic in the treatment of concepts such as interference 
or coherence that are ubiquituous in classical optics. 

However, when one tries to extend this variable into the quantum domain one is immediately faced with two serious 
problems: first, there is no precise prescription to deal with the phase of a three-level atom iQ; second, we lack 
a satisfactory description of the atom-field interaction in terms of relative phase. The main goal of this paper is 
precisely to remedy these drawbacks by proposing a comprehensive theory of the relative phase for three-level atoms 
interacting with quantum fields. 

Phase for three-level systems has been handled by invoking fuzzy concepts such as the phase of the associated 
wavefunction We emphasize that these notions, though well established in the classical limit, are not easily 

extrapolated into the realm of the quantum world. Since phase is considered to be a physical property, in the 
orthodox picture of quantum mechanics it must be associated with a selfadjoint operator or at least with a family of 
positive operator-valued measures p^ . In this spirit, phase operators for the algebra su(2), which describes two-level 
systems, have been previously worked out Since su(3) is the natural arena in which to deal with three- level 

atoms, we propose here an extension enabling us to introduce phase operators for these systems with a clear physical 
picture. 

Concerning the aforementioned second problem, we stress that an operator describing the relative phase between 
field and atoms has resisted a quantum description For this reason, its role has been played by the field 

quadratures p3], which in many respects properly account for phase relations and are free from the difficulties that 
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phase encounters in the quantum domain. When focusing on the relative phase between two subsystems, we think 
the best way to proceed, much in the spirit of our previous work on the subject |2^, is to try a polar decomposition 
of the quantum amplitudes that parallels as much as possible the corresponding classical factorization. 

For the relative phase between two quantum field modes this is quite a straightforward procedure and leads to a 
unitary solution p^ . For the case at hand, this polar decomposition seems to be more involved, mainly because, 
unlike for the case of two harmonic oscillators, the Hamiltonian cannot be cast in terms of su(3) generators, but 
rather in terms of a polynomial deformation of su(3). These nonlinear algebras have been examined very recently in 
quite different physical contexts ]25[| . Without embarking on mathematical subtleties, we shall exploit these previous 
results to perform the polar decomposition in an elegant way, obtaining a bona fide Hermitian operator representing 
the relative phase we wish to examine. 

It is our firm belief that this Hermitian operator and its eigenstates should be a fundamental tool for examining 
all the issues related with quantum interference between atomic pathways. This leads to coherent effects as those 
mentioned earlier. In this paper we use this operator to introduce its associated probability distribution and to explain 
some relevant dynamical features of the model. Nevertheless we stress that given the variety of phenomena of both 
conceptual and practical importance arising in this field, it is outside the scope of this paper to provide a full account 
of them. Rather, our goal is to emphasize the role that the relative phase plays in atom-field interactions, a role that 
has previously gone almost unnoticed. 

II. PHASE FOR THREE-LEVEL SYSTEMS 

We wish to explore in some detail the phase properties of a three-level system. To be specific we shall consider a 
A configuration, as shown in Fig. 1, with energy levels loi < L02 < (we shall use throughout all this paper units 
h — \) and with allowed dipole transitions 1 3 and 2 <-> 3, but not 1 <-> 2. 

In order to describe the kinematics of this system one must take into account that the natural generalization of 
the Bloch vector [ p2| (which is essential in understanding the behavior of a two-level system) comprises now the nine 
operators 

^'^■-b'X^I, (1) 

where |z) denotes the eigenstate of the ith atomic level. One can easily check that they satisfy 

= (2) 

which correspond to the commutation relations of the algebra u(3) |26| . 

Obviously, the three "diagonal" operators iS" measure level populations, while the "off-diagonal" ladder operators 
S'*-' generate transitions from level i to level j. To emphasize this idea we define raising and lowering operators by 

S% = S'^ when j > i, 

(3) 

S'^_l = when j <i. 

Because one has the trivial constraint S"^^ -|- S'^'^ + S^^ — I, only two populations can vary independently. For this 
reason, it is customary to introduce two independent traceless operators 

413 = 1(^33 _ ^11)^ 5f = i(^33_^22)^ (4) 

that measure atomic inversions between the corresponding levels. The operators (5*^?, Sl^) (i 7^ j) turn out to be the 
eight generators of the algebra su(3). Note that {S]?, Sl^) and (S'^3 5*23-) foj-^j two su(2) subalgebras that physically 
correspond to the two dipoles that appear in the two allowed transitions 1^3 and 2^3. However, these two dipoles 
are not independent, since Eq. (^) imposes 

[5i^^23] = -5i^ [Sl^Sf] = 0, (5) 

which constitutes a highly nontrivial coupling between them. 

It is usual to take for granted the existence of the phase of the atomic dipoles, mainly because classically one has a 
clear picture of their meaning. However, strictly speaking, we do not have any prescription to deal with this variable 
at the quantum level. Nevertheless, as anticipated in the Introduction, one could expect to get a formalism very close 
to the one developed recently for two-level systems. 
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For definiteness let us first focus on the transition 1^3. If only these two levels were involved, such a transition 
would be described by a superposition of the form 

I*) = sin(i?/2) |1) + e"^ cos(?9/2) |3). (6) 

This corresponds to a 1/2 angular-momentum system. For this case the Bloch vector is S = o-/2, a being the 
Pauli matrices, and it is easy to work out that the mean values Sj — {^\Sj\^) are given by 

Sx = sini? cos tp, 
Sy = sin I? sin 

= cos??. (7) 

This would support the naive belief that, when viewed in the Bloch sphere s'^ + Sy + — 1, the parameter is the 
phase angle associated with the atomic dipole 1 <-> 3 and is canonically conjugate to Sz p8| . In consequence, there is 
a widespread usage of dealing with this dipole phase ip as a state parameter instead of a quantum variable. To gain 
further insight into this point, let us note that 

S- = Sx — isy = sin-i? e^'^ , (8) 

so it is clear that the parameter "dipole phase" can be obtained through the decomposition of s_ in terms of modulus 
and phase. Obviously, it is tempting to pursue this simple picture by taking into account that, at the operator 
level, the equivalent to the decomposition in terms of modulus and phase is a polar decomposition. Thus, it seems 
appropriate to define the operator counterpart of Eq. dq) in the three-level case we are considering as 



51^ = ^Sl^Sl^ E^^. (9) 

Here E}^ — exp(i(^^'^), where (p^^ is the Hermitian operator representing the phase of the transition. 
One can work out that the unitary solution of Eq. (|^) may be written as 

El^ = \l){i\+e'^'o |3)(l|-e-^'^o |2)(2|, (10) 

where the undefined factor e"^o appears due to the unitarity requirement of E^ . The main features of this operator 
are largely independent of (^q, but for the sake of concreteness, we can make a definite choice. For example | |20| , if we 
assume a state such as the linear superposition (||), the complex conjugation of the wavefunction should reverse the 
sign of 1^^'^, which immediately leads to the condition e^'^o — — i. We conclude then that the operator we are looking 
for can be represented as 



The associated eigenstates are 



^i3 = |l)(3|~|3)(l| + |2)(2|. (11) 
l-^o') = |2), 



(12) 



l^i?) = i=(|3)±z|l)), 



and the eigenvalues of i^^'^ are and ±7r/2, respectively. This is a remarkable result. It shows that the eigenvectors 
|iy9j?) look like the standard ones for a quantum dipole or a spin 1/2. However, the "spectator" level |2) (i.e., it 
apparently does not take part in the transition) is an eigenstate of this operator, which introduces drastic changes in 
the dynamics. In other words, the phase of the transition 1^3 "feels" the state |2). 

Perhaps, the most striking feature of this theory is that a measurement of the phase of the transition gives only 
three possible values: and ±7r/2. While this kind of statement seems rather reasonable when dealing with spin 
systems, they are scarcely recognized when dealing with atoms. 

An analogous reasoning for the transition 2 3 gives the corresponding operator E"^ 

^^3^|2)(3|-|3)(2| + |1)(1|, (13) 
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with eigenvectors 



and the same eigenvalues as before. 
Note that 



(14) 

Wf) = i=(|3)±i|2)), 



El^Ef^^El\ (15) 



which clearly displays the quantum nature of this phase. 

To any smooth function F{if^^) of the dipole phase ip^^ we can associate the operator 

F{^^^)= J2 n^i') iv'r'i (16) 

where the sum runs over the eigenvalues and ±tt/2. The mean value of this operator function can be computed as 

{F{0^^)) =J2f{^1^)P{^1^), (17) 

r 

where P{(p^J) is the probabiHty distribution 

p(^;^HTr[^|^;-'-)(^;^l], (is) 

for any state described by the density matrix q. We shall make use of these results in next Section. 

III. EXPLORING THE ROLE OF THE RELATIVE PHASE 

A. Quantum dynamics of a three-level atom coupled to a two-mode field 

To maintain the discussion as self contained as possible we briefly review the basic aspects of the interaction of 
a three-level atom with a two-mode field. We shall merely outline what will be essential for our later discussion of 
phase properties. We write the Hamiltonian for this system as 

H = Ha + Hf + V, (19) 

where 

Ha = (^jS", 

i 

V = gaiaSl^ + S^J) + 9b{bSf + S^J). (20) 

Here Ha describes the dynamics of the free atom and Hp represents the cavity modes of frequency cua and Wb, with 
annihilation operators a and b, respectively. Finally, in the interaction term V, written in the dipole and rotating- 
wave approximations, we assume that the allowed transition 1 3 couples (quasi)resonantly to the mode a and the 
transition 2^3 couples to the mode b with coupling constants ga and gb that will be taken as real numbers. 

The bare basis for the total system is \i)A \na,nb)F, where \na,nb)F, is the usual two-mode Fock basis. However, 
one can check that the two excitation-number operators 

Na^a^a- S^^ + 1, Nb = b^b- S^^ + 1, (21) 



are conserved quantities. In consequence, we can rewrite the Hamiltonian (19) as 

H^Ho + Hint, (22) 
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where 



Hq = LOaNa + UJbNb + {uJs - CJq - UJb)I, 

(23) 

Hint = — AqS*^^ — AftS*^^ 



The detunings are defined as 



AQ=tJ3i-Wa, Ab = W32-a;{„ (24) 
with It is a straightforward calculation to check that 

[^O,^int]=0. (25) 

Therefore, both the free Hamiltonian Hq and the interaction Hamiltonian Hint are constants of motion. Hq determines 
the total energy stored in the system, which is conserved by the interaction. This allows us to factor out exp(— iifoO 
from the evolution operator and drop it altogether. Thus, the problem can be reduced to study the restriction of Hint 
to each subspace Ti'^"'^'') with fixed values of the pair of excitation numbers {Na, Nt). In each one of these subspaces 
'j^(Na.,Nt) fjnQYe are three basis vectors that can be written as 

\i;na ^ Na - fJ.t,nb = Nb - Vi), (26) 

where the values of fii and Vi are defined as 

(m, M2, Ms) - (0,1,1), (z.1, :/2, 1^3) = (1,0,1). (27) 

Note that when Na = 1 and A^b = or Na = and A^f, = 1 some states may have negative photon occupation number 
and must be eliminated. In the subspace Ti^^'^'^''\ Hint is represented by the 3x3 matrix 



/ gbVNb gaVNa \ 
H^^'r'''^ = gbVm ~Ab . (28) 

\gaVN: -Aa J 

Let us assume that at i = the atomic wave function can be written as the coherent (normalized) superposition 

|vI'(0))a = ^q |z)a, (29) 

i 

while the field is in a two-mode coherent state 

|*(0))f = |aa,a6) = Qn^n, \na,nh)F- (30) 

Here Qn is the Poissonian weighting factor of the coherent state (with zero phase) with mean number of photons n: 



(31) 

At a later time t the state vector for the atom-field system in the interaction picture can be expressed as 

OQ 3 

Na,Nb=0 i,j=l 

X ^f-^^)(0 \i,Na-ti^,Nb~iy,), (32) 
where Z^^^'"^''^(i) are the matrix elements of the evolution operator in the subspace H.i^'^'^b) 

u[^-''''\t) = {i- Na-^l^,Nb- I eM-^Hl^t^'^^'^i] b'; Na - ,Nb-y.j), (33) 

which can be calculated exactly as can be seen, e.g. in Ref [|l2j. This state describes completely the system evolution 
and will be the basis for our phase analysis in the following. 
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B. Polar decomposition of the relative atom-field amplitudes: deformed su(3) dynamics 

As stated in the Introduction, our goal is to describe the atom-field relative phase by resorting to a polar decom- 
position of the corresponding complex amplitudes. To this end, let us define the operators 



yl3 _ s cl3 _ cl3. 



= bSf, xf = sf. 



(34) 



These operators satisfy most of the usual su(3) commutation relations provided ^ is recast as 

j^l3^^23] ^ _^_12^ (35) 

where 

However, some of them must be modified in the following way 

= 7Va(l-2S'"-5^2), 

which corresponds to a polynomial deformation of the algebra su(3). The essential point for us is that one can develop 
a theory in very close analogy with the standard su(3) algebra. 

For simplicity, let us focus first on the allowed transition 1 <-> 3 and notice that in every three-dimensional invariant 
subspace Ti'^^'"^''-' the state |1; iVa, — 1) plays the role of a vacuum state since 

X"|i;7V,,iVb-l) -0. (38) 

In this subspace the operator X},'^ is diagonal and we can work out again a polar decomposition similar to (|^) , namely 



113 = ^113113 ^13^ (39) 



where the operator \l X^X]^ is diagonal and 

(40) 



[E]^,Na] = [E^i,Nb]=^- 

The first equation ensures that the operator E^ — exp(i$^'^), representing the exponential of the relative phase 
between the field and the dipole 1 <-> 3, is unitary. The second one guarantees that we may study its restriction to 
each invariant subspace Tl'-^^'^''\ 

Much in the same way as we did in Sec. II, the operator E^ solution of ( |39| ) can be expressed in 7^'-^°'^''-' as 

Ei^ = \l;Na,Ni,~l){3;Na-l,Nb-l\ 
- \S;Na~l,Nb-l){l;Na,Ni,-l\ 

+ \2;Na-l,Nt,){2;Na-l,Nb\. (41) 
As one would expect, it acts as a ladder- like operator 

El^\3;Na-l,Nb-l) = |l;iV„7Vfc-l), 

(42) 

Ei^\2;Na~l,N,) = \2;Na-l,Nb), 
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with eigenvectors 



^ -L(|3,7Va-l,iVb-l)±i|2;iV,-l,iVb)), 



(43) 



and thus has eigenvectors 

^ |2;7V,-l,iVb), 

^ -L(|3,iV^_i,JVb-l)±i|l;iV<„iVfc-l)), 

while the eigenvalues of are and ±7r/2, respectively. 

Obviously, a similar reasoning for the transition 2 <-> 3 gives the corresponding operator as 

Ef = \2;Na-l,Ni,){3;Na-l,Ni,-l\ 
- \3;Na-l,Nb-l){2;Na~l,Nb\ 

+ \l;Na,Nb-l){l;Na,Ni,~l\, (44) 



(45) 



and the same eigenvalues as before. The states (^^ and ( p5| ) are the basis for our subsequent analysis of the dynamics 
of the relative phase. 



C. Relative-phase distribution function 



For any state, the information one can reap using a measurement of some observable is given by the statistical 
distribution of the measurement outcomes. Once again, let us first focus on the relative phase between the field mode 
a and the dipole transition 1 <-> 3. According to Eq. (|l^) it seems natural to define the probability distribution 
function of a state described by the density matrix g{t) as 



13X/^13|]^ 



(46) 



where the vectors 1$^^) are given in Eq. ( ^3| ) and the subscript r runs the three possible eigenvalues 0, and ±7r/2. This 
expression can be interpreted as a joint probability distribution for the relative phase and the excitation operators 
Na and Nb. From this function, we can derive the distribution for the relative phase as the marginal distribution 



For a general state as in Eq. (63), one has 



which, through direct calculation, gives 



E 



i=i 



P(*±',t) 



E 



(47) 



(48) 



(49) 
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Much in the same way one also gets analogous results for the transition 2^3: 

2 

3 



(50) 



Na,N,,=0 



i=i 



This is our basic and compact result to analyze the evolution of the relative phase. 

We have numerically evaluated this distribution for the three allowed values of the relative phase for the case when 
the atom is initially in the ground state |1) and modes a and b are in a coherent state with a mean number of photons 
fia and n^, respectively. For computational simplicity we have used the rescaled time 



9 at 



(51) 



in all the plots and have assumed that — gi,, which is not a serious restriction. 

In Fig. 2 we have plotted a typical situation of a weak field, in which the number of excitations in the system is 
small, say ~ fif, ~ 1. The pattern shows an almost oscillatory behavior, which can be easily understood if we 
retain only the two first terms in the sums over Na and Nh and use the explicit form of the evolution operator U. 
A relevant and general feature that is apparent from this figure is that the probabilities associated with and <i>*^ 
always oscillate out of phase, a point previously demonstrated for the case of the Jaynes-Cummings model pof . 

Perhaps more interesting is the case of strong-field dynamics, when the number of excitations in the system is 
large and so or n;,, or both, are large. In Fig. 3 we have plotted the relative-phase probabilities for ha = 50 
and (a) hi, = 0.5, and (b) hh = 50 photons, with the atom initially in the ground state |1). When hi, — 0.5, the 
distribution P(<I>J^, t) (which is the probability of finding the atom in the level |2)) is almost negligible, while P($j?, t) 
show collapses and revivals. One may interpret this physically as follows: the transition 1 <-> 3 is so intense due to 
stimulated processes in mode a that there is no population transfer to level |2), which originates a regular oscillation 
of the dipole 1 3 with the corresponding collapses and revivals in the relative phase. The well-known (nearly) 
time-independent behavior in the time windows between collapse and revival is also clear. The probability of finding 
the atom in the level |1), P(<I>o^,t), tends to be 1/2 (except at the revivals), which confirms that the transition 1^3 
is almost saturated. 

When hf, grows, the position of the collapses and revivals changes, according to standard estimates |jl^. When 
ha — hi, = 50, P{^Q'^,t) is centered at 1/4, while P(<I>g^,<) is centered at 1/2, showing that the populations tend to 
be equidistributed because now the transition 2 <-> 3 is almost saturated too. 

Very interesting physical phenomena arise when one considers coherent superpositions of atomic states, because it 
is then possible to cancel absorption or emission under certain conditions, i.e., the atom is effectively transparent to 
the incident field even in the presence of resonant transitions. A semiclassical analysis p^], in which the fields are 
treated as c-numbers and described by the complex Rabi frequencies fiae"'^" and rif,e~*^"^(note that 6a and 9b are 
the 'phases' of the respective fields), easily shows that when the initial atomic state is a superposition of the two lower 
levels of the form 



coherent trapping occurs whenever 



|*(0)), 



1 



V2 



(|1> 



|2)), 



'^a- Ob - = iTI"- 



(52) 



(53) 



In other words, when these conditions are fulfilled the population is trapped in the lower states and there is no 
absorption. 

To corroborate this behavior valid in the strong- field limit, let us note that, although the transition 1 ^ 2 is 
dynamically forbidden, one can still define phase eigenvectors for it: 



= |3,7V,-l,iVfa-l), 



(54) 



1 

7^' 



l*±) - ^{\2-Na~l,Nb)±i\l;Na,Ni,~l)). 
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In Fig. 4 we have plotted the probabilities P($J^,t) when the atom is initially in a trapped state like ( [S^ ) with 
TT-a = nt — 50. We see that P{^l^,t), which is the probability of finding the atom in the upper state, shows the 
remarkable behavior of trapping, except for the presence of very small superimposed oscillations. Note also that 
P($j?,t) also shows the same kind of behavior. 

IV. CONCLUSIONS 

In this paper we have investigated an appropriate operator for the quantum description of the relative phase in the 
interaction of a three-level atom with quantum fields. To this end, we have resorted to a proper polar decomposition of 
the corresponding amplitudes, which has been justified on physical grounds as well as using the theory of polynomial 
deformations of su(3). 

From the phase states obtained in this procedure we have defined a probability distribution function for the relative 
phase and studied its time evolution, showing how the formalism could be applied to understanding more involved 
phenomena. 
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FIG. 1: The energy scheme of a three- level A- type atom interacting with two single- mode quantum fields, coupling the two 
ground states to a common excited atomic state. 



FIG. 2: The probability distribution function for the six allowed values of the relative phase as a function of the rescaled time 
r in the case of a weak field with fia = nt, = 1. 



FIG. 4: The probability distribution function for the allowed values of the relative phase "1> as a function of the rescaled time 
T for a trapped state with fia = fib = 50. 



FIG. 3: The probability distribution function for the six allowed values of the relative phase as a function of the rescaled time 
T in the strong-field limit with: a) fia = 50, fib = 0.5 and b) fia = 50, fib = 50. 




0.00 J ^ ^ 0.00 

0.00 0.27 0.54 0.81 1.08 0.00 0.27 0.54 0.81 1.08 1.35 



(a) 



+ 

-©- 



o 

-©- 



-©- 

Oh 



0.99 
0.66 
0.33 
0.99 
0.66 
0.33 
0.99 
0.66 
0.33 



t 



0.00 



0.00 0.27 0.54 0.81 1.08 1.35 



